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LLR analysis at IfE - program LUNAR 


■ Goes back to mid of 1980ies (FESG, Munich) 

■ 2 main development phases: ~1985-2001 (Munich), since 2006 at IfE 

■ LLR-only analysis 

- ephemeris computation based on initial values from DE421 (Sun, planets, 
largest asteroids) 

- no further planetary data (optical, radar,...) included 

- ephemeris model: EIH equations of motion for all bodies as point masses 

- Earth orientation: IERS conventions 

- lunar orientation: Euler equations integrated together with translation 
including relativistic corrections (geodetic precession, Lense-Thirring) 

- additional forces due to inhomogeneous gravity fields, tides 

■ Combined analysis of solar system data and LLR planned in future 
(project in research unit FOR 1503) 










Program updates at IfE 


■ Data reduction compared with IMCCE, good agreement 

■ Results show some room for improvements updates in ephemeris 
computation needed 

- many modeling parts still from first development phase 

- simplifications (due to computation time reasons and accuracy 
requirements) 

- slight inconsistencies in force model (,interactions') resolved 

- 3 parts: Earth-related, additional gravitational effects, Moon-related 
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Program updates at IfE - Earth related 


■ Earth tides - tidal acceleration 

- former model 

• single lag angle of time delayed tides 

• only Moon as tide generating body 

• effect on lunar translation 

- new model (according to DE430 ephemeris) 

• degree 2 tidal potential (-98% of tidal effect) 

• arbitrary tide generating body possible 

• 5 tidal time delays (2 estimated, 3 fix at DE-values) 

• effect on Moon: 

■ via change in Earth‘s degree 2 potential coefficients on lunar 
translation and rotation 

■ via acceleration on lunar translation 



r 



EM 


ife) 


DFG Forschergruppe 
, Referenzsysteme 


I j l Leibniz 

( Of Z Universitat 

3 0; 4 Hannover 










Program updates at IfE - Earth related 


■ Earth tides - tidal acceleration 
- some results 

• Sun + Moon as tide generating body 

• tides from Jupiter, Venus -> <0.1 mm in r EM in 45 years 

• estimated time delays 

■ diurnal t 21 =575 s (DE430: 636 s) 

■ semi-diurnal t 22 =226 s (DE430: 219 s) 




■ Secular trend in C 20 

- several models tested (linear, quadratic) 

- best result in LLR analysis with model 
from IERS Conventions 2010 with 

C 20 = 2.6xl0 _11 yr _1 
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Program updates at IfE - additional gravitational effects 


■ Interaction of Sun/planets with Moon 

- former model - interaction of 

• planets with point-mass Moon 

• Sun with lunar degree 2 

- new model 

• planets with lunar degree 2 

• Sun with lunar degree 2+3 

- main effect on lunar rotation, e.g. 

• ~ 19 mm on surface from Venus 

• ~ 4 mm on surface from Jupiter 
(from ephemeris with equal initial conditions) 

- in analysis: residuals decrease ~0.1 mm on average 
-> small effect but maybe needed in future 
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Program updates at IfE - additional gravitational effects 


Figure-figure interaction between Earth and Moon 

- former model 

• simplified version of degree 2 - degree 2 coupling 

• effect on lunar rotation considered 

- new model 

• coupling up to any degree/order of the 
gravitational field of Earth and Moon possible 
(Ilk, 1983) 

• effect on translation and rotation 

- results 

• improvement due to complete 
degree 2-degree 2 coupling 

• Earth degree 2 - Moon degree 3 
-> some mm on surface 
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Program updates at IfE - Moon related 


■ Rotation of deformable Moon 

- former model 
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• tidally deformed tensor of inertia only in rotation 

• no changes of potential coefficients (interaction with Earth, additional 

effect on translation and rotation) ♦ 

• no consistent core-implemention 

- new model (according to 2-layer model of DE430) 

• basis tide-free tensor of ineria, elements from 

■ C 20 (GRAIL) 

■ C 22 , dynamical p (estimated in analysis) 

■ dynamical y derived 

• core moments with DE-fixed values for inertia-ratio C core /C Moon 
and core flattening 
















Program updates at IfE - Moon related 



• tidal + spin deformations on mantle tensor of inertia 

■ Love number k 2 fixed on GRAIL-value 

■ tidal deformation due to Earth (much larger than from Sun) 

■ 1 time delay used 

-> equations for rotation complete 
-> coupled differential equations for core + mantle 
-> coupling constant, initial rotation vector of core estimated 
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• degree-2 changes in potential coefficients from deformed mantle tensor 
-> enters in computation of external forces (translation/rotation) 
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Program updates at IfE - Moon related 


■ Effect of new rotation¬ 
modelling on post-fit 
residuals, wrms (without 
core) 


■ Effect including core 
(which leads to improved 
libration modelling) on 

2-way residuals 
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Results - wrms 


■ 1 -way wrms, comparison with former model 



■ Further investigations 

- effects in longitude libration visible ^empirical correction as in DE430 

- core flattening not yet estimated -> possible effect on librations 

- tidal and lunar rotation modelling -> room for improvement 
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Results - coefficients lunisolar nutation 


i J v \ u-- / 


Periode MHB2000 [mas] old model new model (two test cases) 




less accurate, leads 
to larger discrepancies 
in other coefficients 
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Results - relativistic parameters 


■ Temporal variation of gravitational constant 


- modelled as G(t) = G 0 
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- moderate to strong correlations with 

• lunar core rotation vector (fixed in solution) 

• some station coordinates (constrained a little bit) 

- as single parameters: — = (0.7±0.8)xl0" 13 yr _1 

G 0 
• • 

— = (1.6±2.0)xl0" 15 yr" 2 

G 0 

• • 

c G 

- estimated together: — = (0.8±l.l)xio _13 yr 1 — = (-0.3±2.4)xl0" 15 yr" 

G 0 G 0 
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Results - relativistic parameters 


■ Equivalence principle 

- estimating ratio A(m /m,) 

V 9 '/EM 

• partials: computed numerically by introducing 
additional acceleration of Moon towards Sun 

• Afm/rn) = (-3.0 + 6.6)xl0“ 14 

V 9 1 /EM 

- Estimating Nordtvedt parameter rj 

• partials: analytical from synodic range 
variation: 13.1m*cos(D)* rj 

• r| = (-0.2+1.2) x 1CT 4 



- into direction of galactic center (e.g. due to dark matter) 

• amplitude: A =3.0±3.3mm 

9 C 

• additional acceleration a =(-l.l±1.2)xl0" 6 
in parts of 1.9e-8 cm/s 2 9C 
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Results - relativistic parameters 


■ PPN parameters p, y 

- included in EIH-equations of motion, partials numerically 

- correlated with station coordinates (constrained) 

- show also correlations with z-coordinate of lunar initial orbit values 

- p—1 = (0.9±1.0)xl0 -4 

- y _] = (_i_ 2 ± 1.6)xio -4 (not as accurate as Cassini-result) 

- p from combination of PPN-parameters 0.25(y+r|+3) and Cassini- y 
p-l = (0.03±6.1)xicr 5 

■ PPN preferred frame a,, a 2 w.r.t. cosmic microwave background 

- a, =(-l.l±2.0)xl0‘ 5 

■ a, =(-0.6+0.9)xl(T 5 ( n °f as accurate as test with Sun's spin) 










Results - relativistic parameters 


■ Geodetic precession of lunar orbit 

- introducing GP a second time as additional acceleration 

- factor h gives relative deviation in from Einstein‘s theory (~1.9 as/cy) 

- strong correlation with 

• lunar core rotation vector (fixed) 

• z-component of lunar initial velocity (fixed) 

h = (-0.6 ± 2.0) x 10” 3 

■ Yukawa-term (1/r 2 -test), acceleration due to v^ M — — GMeMm 



interacting range X = 380000km 
coupling constant a estimated 
correlations and fixed values like GP 

a = (-4.0 + 5.0) x 10 -12 











Summary 


■ IfE-LLR ephemeris model updated 

- tidal acceleration, secular trend in Earth s C 20 

- additional gravitational interactions planets-Moon, Earth-Moon 

- lunar rotation as 2-layer core/mantle model 

- effect of lunar deformation in all lunar equations of motion (translation 
plus rotation) 

■ Smaller residuals and more accurate parameter estimation 

- increased accuracy in relativistic parameters (strong limits for validity of 
equivalence principle and gravitational constant) 

- no deviation from Einstein‘s theory of gravity up to now 









